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Abstract 

<^ ■ We compute next-to-next-to-leading order spin contributions to the post-Newtonian equations 

. of motion for binaries of compact objects, such as black holes or neutron stars. For maximally 

O \ spinning black holes, those contributions are of third-and-a-half post-Newtonian (3.5PN) order, 

■ improving our knowledge of the equations of motion, already known for non-spinning objects up 
l/^ ■ to this order. Building on previous work, we represent the rotation of the two bodies using a 

pole-dipole matter stress-energy tensor, and iterate Einstein's field equations for a set of potentials 
parametrizing the metric in harmonic coordinates. Checks of the result include the existence of a 
conserved energy, the approximate global Lorentz invariance of the equations of motion in harmonic 
^ ' coordinates, and the recovery of the motion of a spinning object on a Kerr background in the test- 

mass limit. We verified the existence of a contact transformation, together with a redefinition of 
the spin variables that makes our result equivalent to a previously published reduced Hamiltonian, 
. obtained from the Arnowitt-Deser-Misner (ADM) formalism. 
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I. INTRODUCTION 



Gravitational wave detectors will soon enter a new era, with the advanced versions of 
Virgo and LIGO expected to start operating around 2015, and the construction of KAGRA 
[l|-|3|. Primary targets for these laser ground-based interferometers, or for possible future 
space-based interferometers, are inspiralling binaries of compact objects — neutron stars 
and/or black holes — , which emit gravitational radiation during their late inspiral, merger, 
and ringdown phases. Hunting for the faint signal and best separating it from the noise using 
matched-filtering techniques {e.g. 0, Isf) requires a very precise modelling of the expected 
waveform, which the post-Newtonian (hereafter PN) approximation aims at providing for 
the inspiral phase jsj. 

Up to now, such high-precision PN templates are available for the non-spinning case up 
to the 3.5PN order, i.e. the order in the formal expansion in powers of (with c 
being the speed of light). However, observational evidence points toward the existence of 
fast-rotating black holes, endowed with significant an gula r momentum (or spin), for stellar 
size j7-ll| as well as for supermassive black holes 12-14 . 



It is thus important to complete our knowledge of the PN predictions for the gravitational 
waves emitted by such systems by including the spins of each of the two bodies. This implies 
first taking those spins into account in the dynamics of the binary. For maximally spinning 
objects, the leading-order linear- in-spin (which we shall call spin-orbit) effect arises at 1.5PN 
order when regarding the spin as a 0.5PN quantity [see Eq. (11.11) below]. The next-to-leading 
and next-to-next-to-leading contributions show up at respectively 2.5PN and 3.5PN orders. 
The effect of spacetime curvature on the motion of a spinning test particle was obtained in the 
seminal work by Papapetrou 15[4l7l | , after an earlier derivation by Mathisson (republished 



in [18[). Barker and O'Connell [19|, |20| determined the leading-order spin-orbit and spin- 
spin effects in the two-body dynamics. More recently, Kidder, Will and Wiseman [2l|, 2^ 
computed the corresponding contributions to the radiation field and, most importantly, to 
the orbital phase of the binary, to which the templates are crucially sensitive. Effective field 
theory methods [1^ were also used to rederive the leading-order spin-orbit and spin-spin 
contributions to the dynamics (13] . The problem in the limit _of_a spinning test particle on 

(see also 



a Kerr background has been addressed in 
neglecting the gravitation damping force). 



[251,1261 



27 for a Hamiltonian model 



By means of a PN iteration of the Einstein field equations in harmonic coordinates, the 
equation of motion at the next-to-leading order was first investigated by Tagoshi, Ohashi and 
Owen L 



29| . It was then confirmed and completed (as well as extended to the radiation 



field) by Blanchet, Buonanno and Faye [30|, l3l|- The results for the evolution equations 
were retrieved by two independent calculations, using a Hamiltonian approach in ADM 
coordinates [12] on the one hand, and using effective field theory methods [s^, 34] on the 
other hand. The ADM computation was later generalized to the many-body problem j35| and 
extended to the next-to- leading order spinl-spin2 and spin square {e.g. spinl-spinl) effects 
in Refs. L 
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Next-to-next-to-leading order spin-orbit effects in the Hamiltonian of the 
binary were first computed in [3^ , resorting to the ADM scheme adapted to matter sources 
composed of spinning point particles [11]. In a subsequent work, spinl-spin2 interactions 
terms were also added at the 4PN order [i^. Effective field theory methods progressed 
concurrently by computing the 3PN spinl-spin2 and spinl-spinl contributions 4l|-|43j . and 
the 4PN spinl-spin2 interactions 44 . 
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Our aim in this work is to extend the approach of Ref . j30| jl| based on the PN expansion 
of the metric in harmonic coordinates, in order to compute the spin-orbit next-to-next-to- 
leading order 3.5PN contributions to the equations of motion. We shall thus provide an 



essential validation of the result obtained by the ADM method [38|. Our final purpose is to 



derive the spin-orbit terms in the emitted GW ener gy f lux and in the phase of the binary at 
the same approximation, extending therefore Ref. |3l| to next-to-next-to-leading order, so 
as to provide more accurate PN predictions to the template-based data analysis. 

We adopt the same convention for the PN order counting of the spins as in Paper I; 
namely, we redefine our spin variable with respect to the "true" spin angular momentum 
following : 

S = cS true = Gm^odyX , (1-1) 

where mbody is the mass and x the dimensionless spin parameter (G denotes Newton's 
constant). For a maximally rotating compact object, we have x ~ 1 and S ~ Gm^^^y, so 
that our spin variable can be counted as of Newtonian order in this case. For slowly rotating 
objects however, we have x ~ "^surf/c where Wsurf is the rotation velocity of the body surface, 
which implies that S acquires an additional factor 1/c. In that case, spin-orbit contributions 
effectively appear at 2PN order, and spin-spin ones at 3PN only. Our counting will always 
assume rapid rotation. The leading-order spin-orbit terms will thus carry an explicit factor 
in the present article we shall obtain the next-to-leading and next-to-next-to-leading 
spin-orbit corrections ~ and ~ 

The paper is organized as follows. In Sec. Ull we review the pole-dipole formalism for mod- 
elling compact objects as point particles with spins. In Sec. Illlt we describe the parametriza- 
tion of the PN metric in terms of a set of elementary potentials and present the equations 
of motion as deduced generically from metric components. The core of our work consists 
of the computation of the required potentials in Sees. [IV] and El first using Hadamard's 
regularization (as well as the "pure Hadamard- Schwartz" prescription [iH), and next com- 
pleting some of our calculations by means of the more powerful dimensional regularization. 
Sec. |Vl] contains our results for the precession equations and the acceleration. In Sec. IVIIt 
we explain the various checks that have been performed to validate them: construction of a 
conserved energy, check of the global Lorentz invariance, test-mass limit, and recovery of the 



results of Ref. 38| in the ADM Hamiltonian formulation. We conclude shortly in Sec. IVIIII 



II. EFFECTIVE POLE-DIPOLE FORMALISM 



The starting point of our calculations is the model of pole-dipole particles developed b 



5^ 



Mathisson [isj, Papapetrou 16, 17|, Tulczyjew 46, 43| and generalized by Dixon |48 
and Bailey & Israel |51|. This model allows an effective description of the dynamics of 
bodies that accounts for their spin angular momentum by means of singular Dirac delta- 
functions, making analytical computations tractable. We use the same version of the model 
as in Paper I, staying linear in the spins, but we summarize here the main formulae for 
completeness and refer to Paper I for more details. 

The central assumption is that each particle is described by a stress-energy tensor made 
of two parts, a monopolar one and a dipolar one: T'^'^ = + T^^ . These are built with 
respectively a Dirac delta function and a gradient of a delta function, integrated over the 



^ Wc will refer to Ref. [30| as Paper I throughout this work. 
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world line of the particle, according to: 



TiT ^ r ,r p^^u^^'-^^lpp^ , (2,la) 



oo 



r> + 00 

rr = -c I dr V. 



(2.1b) 



Here r is the proper time measured along the world line, described itself by the particle 
position y^{T); S^^^ denotes the four-dimensional Dirac delta function, g stands for the 
determinant of the spacetime metric gfj,^] = d?/^/(cdr) is the four- velocity of the particle 
(satisfying u^u'^ = —1), its linear momentum, and 5''^'' an antisymmetric tensor that 
represents the spin of the particle. 

The linear momentum and spin tensor S'^'^ obey the Mathisson-Papapetrou equations 
of evolution: 



dr 



'■{p'^u" -p'u'') , (2.2a) 



Do'' 1 

= --R^^^XS''' , (2.2b) 

with D/dr = cu'^'Vu, and R^^p^j denoting the Riemann tensor. As is well known, a choice of 
supplementary spin condition (thereafter SSC) is necessary, in order to obtain the correct 
number of degrees of freedom (see js^] for a summary of the various choices in use in the 
litterature). We adopt the following covariant SSC: 

S^^'p, = . (2.3) 

With this relation in hand, one can combine the equations (12. 2p to obtain the link between 
the four- velocity and the linear momentum p^\ hence we deduce the conservation laws: 

Dm , , 

for the mass and for the magnitude of the spin defined by rn^c^ = —p^p^ and 5*^ = S^'^S^^v/'^ 
respectively, and thus conserved along the particle's world line. 

We now restrict the evolution equations themselves to linear order in spins, neglecting 
any term (9(5'^), quadratic or of higher order. One can readily check the proportionality 
relation between and p^ in that approximation: 

pt" = mcu'' + 0{S^) . (2.5) 

In particular, the SSC now reads S^'^Uy = (9(5*^). Notice that the monopolar part of the 
stress-energy tensor reduces to the one of an ordinary point-like particle, while in general it 
could depend on the spin as well. Equations (12. 2p become: 

D.9'''' 

= OiS') , (2.6a) 
mc—- = --R^n'^S'"' + 0{S^) , (2.6b) 



dr 2 
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We see that in the hnear approximation, the spin tensor is parallely transported along the 
motion, which is actually non-geodesic due to the coupling of the spin tensor with the 
Riemannian curvature. 

In the following, we shall use the 3-dimensional form of the energy-momentum ten- 
sor f l2.ip . Using a 3 + 1 splitting of spacetime, the particle's position and coordinate 
velocity are denoted y'^ = (ct,y(t)) and v^{t) = (c, v(t)) (where = cu^/vP, with 
vP = 1 / \/ —gpaVV / c^) , and the spin tensor S^^{t) is considered a function of time. From 
now on, we use boldface letters to denote three-dimensional vectors. We have 



, ^(^)(x-y(t)) ^ 



(2.7a) 
(2.7b) 



where 5^^^ is the three-dimensional Dirac delta function. Expliciting the covariant derivative 
in the dipolar stress-energy tensor, we obtain an alternative form in terms of ordinary 
Christoffel symbols: 



-9TI; 



-a 



SP(^'v^'^S^^)(x - y(t))l + Sf^''T';lv''S^^\x - y(t)) 



(2.8) 



Note that in f l2.7p -( l2l8i) the factor has the generic field point (ct,x) for argument. 

Because of the SSC, we may work only with the spatial components S^^ of the spin tensor, 
and eliminate the 5"°* components according to: 



(2.9) 



In Paper I, a spin vector was used instead of a spin tensor. First, a spin covector was 
defined, consistently with the SSC S'^'Uy = 0{S^), by: 

S^"" = e'^'P'^UpS, + C(5^) , (2.10) 

where e^'^'"^ is the Levi-Civita tensor such that e^^"^^ = —1/ \/—9- Then, a spatial spin vector 
S'pBB constructed as S^^-^ = Y''Sj, with 7*-' being the inverse of the spatial part of 
the metric, i.e. Y''9kj = ^^j- In the present paper we shall mostly work with the spatial 
components S^^ of the spin tensor. This presents two advantages: first, it somewhat simplifies 
the algebra and the index structure by getting rid of the Levi-Civita tensors, and second, 
it makes the check of the Lorentz invariance more straightforward, since S^^ is directly the 
spatial components of the tensor S'^'^. We provide in Appendix [X] the explicit link between 
S''^ and the spin variable Sp^^. 



III. POST-NEWTONIAN METRIC AND EQUATIONS OF MOTION 

As in Paper I, we use harmonic (or De Bonder) coordinates, defined by the gauge con- 
dition du{\/—99'^'') = 0. The PN iteration of Einstein's field equations developed up to 
3.5PN order in Ref. (ssf is valid for a generic matter source, the only hypothesis being that 
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the stress-energy tensor T'^^ must have compact support. We will need here the full 3.5PN 
metric, which is parametrized in terms of elementary potentials asjl 



:V 



X + ViV, + 



6 



32 



+ — T - -VX + R,Vi - -VVy, - 777 + 



9oi 



-Ri 



48 



0{9) 



(3.1a) 
(3.1b) 



9ij = 



1 + + ^V' + - { X + V,V, + ^ 

- 0{8) . 



16 



+ -Wii + — ( + -m, - ViV, 



(3.1c) 



From now on, spatial indices i,j, . . . will be raised and lowered using the Kronecker metric 
6ij, and we will write them indifferently as upper or lower indices. The potentials are defined 
by means of the usual retarded (7^) inverse flat d'Alembertian operator. 



(□7^V)(x,t) 



An 



d^x' 



X — X' 



7/ x',t 



X — X 



where the sources / are built with the following matter quantities. 



cr 



(Tj = -I 

C 



(3.2) 



(3.3) 



as well as with products of derivatives of lower-order potentials. The potentials required for 
the 2.5PN equations of motion were already used in Paper I and read: 



V = □^i[-47rGcT] , 



X = D^' 



AirG Van + WiAiV + 2VidtdiV + Vd^V 



Ri — 



-AttG {Vai - Via) - 2dkVdiVk - -dtVdiV 



= [-AnG [a,, - 5,jauk) - d,Vd,V] , 
while the potentials required for the 3.5PN order are given the following definitions: 



(3.4a) 
(3.4b) 



(3.4c) 

(3.4d) 
(3.4e) 



-A.G[^.,W,^\v^an^aV^.]^Z.A.V^RAd^ 



2d,VjdjRi - d,V,dtW,, + Wdtd.V + 2VidjVidjV + -VAVd.V 



We denote by 0{n) remainder terms of order (n/2)-PN. i.e. behaving as 0(l/c"). 
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(3.4f) 



1.1 1 
-A-kG { -aRi - aVVi + -akWik + -aikVk + -akkVi ] + WkidkiVi 



2 



-':;,vdtVd,v - 2VdkVdkV, + vd^v, + 2VudkdtV, 



-2d{iVkdkVj) - 6ijdkVm{dkVm - dmVk) - -^5ij{dtVf 



(3.4g) 



(3.4h) 



Notice the difference of structure between the sources of these potentials: some sources are 
proportional to one of the a quantities and are compact-supported, while others are only 
proportional to metric potentials. We will call the latter non-compact supported, since their 
source extends in all space. 

Next, in keeping with notations of Paper I, we rewrite the covariant equation of motion 
(I2.6bp in the following 3 + 1 form: 



-TT^ = F^ + ^^ 

at 



(3.5) 



The leading order of the left-hand side of the force law (13.51) is simply the ordinary acceler- 
ation a* = d^ydi. The Newtonian-like linear momentum and forces are 



Pi = giuU^v 



2mc 



(3.6a) 
(3.6b) 

(3.6c) 



With transparent meaning we shall often call Fi the "geodesic part" or the force law (or 
rather the acceleration), while J-^ will be referred to as the "Papapetrou part", deviating 
from geodesic motion. 

Expanding the metric in terms of the elementary potentials (13.41) . one gets the following 
3.5PN expressions for Pi and Fi, which would correspond to the coordinate acceleration of 
a particle following a geodesic motion: 



+ ^ ( + 3Vv' - AVi 



- AVVi + AWijV^ - 8R, 

2 



1 ( + ^ V^; V - - GV.VVV^ + -V^ y V 

c6 Vl6 8 2 ^ 4 

+ 12Xv' + IQZijV^ - lOV^Vi 



-SWijVj - 8VRi - IQYij + C(8) , (3.7a) 
= d^V 

+ -i [i^i^ - '^^i^j ^^^^ + ^VdiVv^ + 2diWjk v^v^ - Wj diV 



- W diVj - SdiRj + ^V^ diV + 8Vj d,Vj + AdiX 



1/11 3 4Q 

\16 2 8 

27 

-lOV,- diV v'^tP - IWdiVj v^v^ - idiRj v^v^ + —V^ diV 
+12VAVi v'^ + QW.k diV v^v^ + W diW.u v^v'' + eaXt;^ 



+SdiZ.k vh^ - 20ViVdiVv^ - lOV^ diV^ - 8Vk diW.k 



Ji^jk v 'iJ — VjVdi V V — iu V Cj \ j V — o Vfc Cj rr jfc ' 
-%Wik diVk - 8Rj diV - 81/ diRj - IQdiYj 

-^V^ diV - 4Vj Vj diV + 16^^ diVj + IQVj diRj 



-WVj diVj - AXdiV - AVdiX + 16^^T j + C(8) . (3.7b) 

Specializing now to our problem, when computing for instance the equations of motion 
for the body 1, the velocity is to be replaced by v{ and the right-hand sides are to 
be evaluated at yi. However, beware that the metric potentials are generically singular 
at the location of the two particles 1 and 2 where they are meant to be evaluated. This 
evaluation is thus give n a sense through the Hadamard "partie finie" regularization procedure 
explained in Ref. [HI, which is of course to be performed after computing the derivatives 
of the potentials. We adopt the so-called pure Hadamard-Schwartz prescription jisj for the 
practical implementation of this regularization. In particular, we use the distributive rule 
for computing regularizations of products of potentials in Eqs. f l3.7p . writing 

(AB), = {A),{B), , (3.8) 

for A, B among (possibly derivatives of) metric potentials. The other ingredient of this 
regularization is the Schwartz distributional derivative which is used for source terms to be 
integrated in Eq. f l3.2p . in the form of the Gel'Fand-Shilov formula valid for homogeneous 
functions, as given by Eqs. f l4.3p - fl4.4p below. 

In this work, we restrict ourselves to contributions to the coordinate acceleration that 
are of linear order in spins. Considering Eq. (13. 5p . these linear-in-spin contributions to the 
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geodesic part of the acceleration can have two origins at this leveh first, from the spin 
contributions to the elementary metric potentials (13. 4p . which we will have to compute at 
the required order, and secondly, from the time derivative dPi/dt, since it is understood that 
order-by-order replacements of the accelerations are to be performed, which include in turn 
spin contributions starting at 0{3). Table |T] indicates which spin parts of metric potentials 
are needed at which order, for the present computation of the next-to-next-to-leading order 
spin-orbit contributions to the equations of motion. As explained in Sec. HV] some potentials 
were computed in all space, but for the most non-linear ones, only their regularized values 
at the location of one of the particles could be computed. In the latter case, the calculation 
is different for each different derivative structure of the potential, and Table |T] gives the list 
of these needed derivatives. 

On the other hand, the effective Papapetrou part of the acceleration or force J-i, which 
corresponds to non-geodesic motion, can be found in Paper I up to 2.5PN order in terms of 
the spin variable •S'pgg. We give here its complete expression up to 3.5PN order, in terms of 
the spin tensor S^^ and the various metric potentials. Defining 

= ^/s + + ^/r + 0{9) , (3.9) 

and expanding the Riemann tensor in Eq. (I3.6cp in terms of metric components, themselves 
expanded in terms of the elementary potentials (13. 4p . we get: 

fi = S'^ {dtd.V + v'^d.kV) + S^'' {2v^d,,V - 2d,,Vk) , (3.10a) 
fi = S'^ {-v'^djVdkV + 2Vv''dkdjV + 2djVkdkV + 2v^dkVdkV - 2dkVjdkV 
+djVdtV + 2VdtdjV + v^v'^dtdkV + v^d'^V) 

+ S^' (Av^d.VdkV - AdjVdkV, - Ad.VdkV, + Ad.kRj - WvW.kV 

+ AV.dikV + 2v^v'dikVi - 2v'dikWi, - 2v^v'duV, + 2v'dumj 

-2v^mVk - 2v^dtdkV, + 2dtdkmj) , (3.10b) 
= S'^ (AdjRkdkV - 2Vv^djVdkV + WkdjVdkV + Sv'^djVidkVi + 2V\^''dkjV 

+ Sv^'VidkjVi + Av^'dkjX - AdkVdkRj + AVv^dkVdkV - 4V,dkVd,V 

- 2v'djW,idiV + 2iPv'dkVidiV - 2v^dkWjidiV - 2v^v'diVkdiV + 2v''diWkjdiV 

TABLE I. First line: the different spin parts (S) of potentials to be computed. Second line: the 
highest PN order required for the spin part of the potentials. Third line: for which potentials an 
all-space (A.S.) computation of the non-compact support part is possible with known techniques; 
and, if not, which are the derivatives of potentials that have to be computed directly using the 
regularization in 1. In the latter case we employ the technique of regularized Poisson integrals 
described in Sec. IIVDI Note that all-space computations of and Rf at the previous PN order 
0(1) are possible and have been checked to be consistent with the regularized versions. 



Potential 




vf 


1^ 






fs 






Order 


0{7) 




0(3) 


0(3) 


0(3) 


0{l) 


0{l) 


0(1) 


Computation 


A.S. 


A.S. 


{^^X)^ 


{Ri)i, {djRi)i 


A.S. 


(5.T)i 




A.S. 
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+ 2VdjVdtV - 2v''djVkdtV + vh'^dkVdtV - 2v^dkVjdtV + {dtVf 
+ 9,djVkdtVk + Av^duVdtVu - 2dkVdtWjk + 2V^dtdjV + SVkdtdjVk 
+AdtdjX + 6Vv^v''dtdkV - Av%dtdkV + Wv^d'^tV - 4V,-5,V) 

+ S^'' [-SVdiVkdjV - SdiVdkRj + SVv^diVdkV - SVjd.VdkV - Av^v^diVidkV 

+ SdjRAV + SVd^VdkVj + Sv'diVidkVj + Av'diVjdkVi - UvWiVidkVi 

- Av'djVidkVi + AdjVidkWii + AdiVidkWji + ^v^v'^u^Rl + mj^k^V 

- AV^v^du^V + WV.dkiV + Av^v'Vidk^V - UWi.dk^V + AV^dk^Vj 

+ Sv%dk,Vj + Wv^v'du^Vi - \^v%d^yi + m.id^yi - Wv'd^Wi^ 

+ Av^v'dkVdiV, - UduV.diV, + Av^v'd.vdiVk - Av'dy.diV^ - Av^dyAv^ 
+ Av'd.vm^ - Ad.VidiW,^ - Avh'dikk, + Av'w,,dikV - Wvh'dikV, 

- 8v%dikVj + AVv'dikmj + 8v'dikZij + Av^dMidiV + Av'dkWudiV 

- Av^diWikdiV - Av^dkVidiVi + Av^diVkdiVi - AdiWkAV, - AdkWudiVj 
+ AdiW^kdiVj + Av^dkVdtVi - 8dkV,dtVi + Av^d.VdtVk + 8d,VAVk 

- AdjVdtW^k - Av'dtdA + Av^VkmV - SVv^mVk + 8V,mVk 
+ 2v^v'dtdiWik - AiPdtdkRi + AWijdtdkV - SVv^dtdkV^ - SVAdkVj 

+AVdtdkWij + 8dtdkZij-2v^v'dtdimk-2v^d^Wik) • (3.10c) 



If we are computing the Papapetrou acceleration of body 1, J-"}, we have to replace in the 
latter expressions S^^ by 5}"' and by v\. Because of the explicit spin factor, at linear order 
in spins only the non-spin parts of the potentials are needed, and 5*2"' does not appear in 
this part of the acceleration. The factor m = mi in Eq. fl3.9p is at the origin of all terms 
depending on the spins by unit mass, i.e. through /mi or 82/^2 in the final results. 
Note that terms with two velocities appear in Eqs. fl3.10p because of the replacement of 5*1 2 
by S'i"'2 according to ( 12. 9p . Again all products of singular potentials have to be regularized 
following the rule of the pure Hadamard-Schwartz regularization, notably the distributivity 
rule f lXH]) FI However we shall find in Sec. |V]that in order to compute the value at 1 of one 
particular potential, namely djkYi, which is especially singular, a priori requires the use of 
dimensional regularization. 

As we said, since the Papapetrou part of the acceleration already includes an explicit 
spin factor, only the non-spin parts of the potentials are required there, and most of these 



have already been computed in previous works such as [53|, |55|. Notable exceptions are the 



second spatial derivatives of the potentials Ri and Yi, which are defined in Eqs. (13. 4p : we 
address the computation of these new regularized potentials in Sees. |IV] and |Vl 



^ Note that wc checked that, at this order, it makes no difference in the final resuh to evaluate the product 
of derivatives of potentials before or after taking the partie finie at 1. 
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IV. HADAMARD REGULARIZATION OF POTENTIALS 



A. Notation and general points 

Throughout the rest of this paper, we use the following notation. The trajectories of the 
two bodies, with masses nii and m2 and spin tensors and S2 , are denoted yi(t) and y2(i), 
and their coordinate velocities vi(t) and V2(t). The inter-body distance is ?"i2 = |yi — y2|5 
and we set n\2 = {vl — Z/2)/|yi ~ y2|- For a generic field point x, we pose ri = |x — yi|, 
n\ = (x* — ?/^)/|x — yi| and similarly for 2. We denote by di, d] and df the partial derivative 
with respect to x*, y\ and y\. The symbol 1 2 means the expression obtained by the 
exchange of the two particles. Lengthy calculations are performed with the help of the 
scientific software Mathematica®, supplemented by the package xTensor 5^ dedicated to 
tensor calculus. 

We systematically expand the retardations inside the retarded integral ( 13. 2p . truncated 
to the appropriate order. For instance, truncating at 0(3), we have: 

-47rn^V(x,t) = j ^^r^/(^''^)-^/ dV9,/(x',t) + ^ j dV|x-x'|9Mx',t) 

+ C(3). (4.1) 

In this expression, time derivatives may be pulled out of the integrals, and it is understood 
that accelerations and time derivatives of the spin are to be replaced by the already computed 
lower order equations of motion and spin precession equations. Note that the second term 
in this expansion has no dependence to the field point x. 

Next, notice that we are using Dirac delta functions in the stress-energy momentum 
tensor (12.71) . while it makes no sense as a distribution a la Schwartz when acting on the 
class of general functions in the problem, which are generically singular at the locations of 
the two particles. We refer to |53|] for a complete discussion of this issue, which was dealt 
with by defining a class of "pseudo-functions" and using Hadamard regularization to define 
the value of a singular function at one of its singular points. In our case, we will give a sense 
to Dirac delta functions inside integrals, following the rule 

d3xF(x)(5i = (F)i , (4.2a) 



where 61 = 5(x — yi) and {F)i means the Hadamard "partie finie" of F at the point yi [54 . 
This rule is extended in an obvious way to derivatives of delta functions, for instance 

j d'^F{^)dA = -{^^F)l . (4.2b) 

Other issues are the problem of distributivity of the partie finie already mentionned, and 
the treatment of the distributional parts of derivatives. Here we shall not use the "extended 
Hadamard regularization" introduced in Ref. [s^]; instead we will apply primarily the "pure 



Hadamard-Schwartz" prescription described in Ref. [45| which is sufficient for most of our 
computations. The pure Hadamard-Schwartz regularization constitutes the core of the most 
powerful and fundamental regularization procedure which is dimensional regularization and 



has been successfully applied to the problem of the equations of motion in Refs. |45|, |57 . 
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Thus, in principle, the result of the pure Hadamard-Schwartz regularization of Poisson- 
type integrals is to be suplemented by a contribution from dimensional regularization which 
appears when the integral develops a pole cx 3) in the dimension. The pole corresponds 

to the appearance of logarithmic divergences in the Hadamard regularization (see Sec. lIVDjl . 
Here, we stick to the pure Hadamard-Schwartz prescription for the computation of the spin 
parts of potentials, where no problematic logarithms appear and where we expect that there 
is no difference between the extended Hadamard or pure Hadamard-Schwartz regularizations 
and the dimensional regularization (as it was the case up to 2.5PN order in the non-spin 
equations of motion |55|). On the other hand, as already mentioned we do use dimensional 
regularization for one particular potential in the Papapetrou part of the equations of motion, 
as discussed in Sec. IVl 

For the distributional part of derivatives of singular functions, we employ a particular 
form of the Gel'Fand-Shilov formula valid for homogeneous functions. Denoting by di the 
full derivative, including the distributional part, and by d°^'^ the ordinary part, which also 
acts on the Dirac delta functions, we get, in the relevant cases of double derivatives and a 
simple function of type jr^ (which is homogeneous of degree — m), 

^^,f = dff + A [dff] + dT'^D, [/] , (4.3) 
where the distributional part is given when £ + m is an odd integer by 

(-r2-(£ + l)!(tel)! [-^] A^^^g(,,„,p^,x+2P-M) 

{t + m)\ 22p(p-l)!(m-2p)!(^+p)! ' ^ ' ^ 

and is zero when ^ + m is evenQ In the present context we apply the formulae fl4.3l) - fl4.4p 
to the expansion of a singular function / around the two singularities yi or y2. 




B. Leading-order potentials 

We first deal with the lowest-order potentials, namely the computation of the leading 
spin-orbit part of the potentials V, Vi and Wij. We give all the results for completeness, 
translating results from Paper I in terms of the spin tensor 5''-' . By expliciting the expression 
of the dipolar part (12.81) of the stress-energy tensor, one may check that the leading order 
for the spin part of the sources a, Ci, and aij is: 

= ^S'^v\dj5i + 1^2 + C(5) , (4.5a) 

af = ^S'^d,5^ + 1^2 + 0(3) , (4.5b) 

a^^ = -ls\^'v{^dA + 1^2 + 0(3) , (4.5c) 

c 

where {ij) indicates symmetrization, 5i stands for 5(x — yi), and the gradients are taken 
with respect to the field point x. Given these expressions for the sources of potentials, and 

Here we pose pq = Max[l, {ra — I ~ l)/2], and [to/2] means the integer part. Notation for multi-indices 
is for instance L = jij2 ■ ■ ■ je, i.e. I is the number of indices on — n^^ ■ ■ ■ n-". We denote = 
^jih ' ■ ■ ^j'>k-ij2k where 2fc = + 2p — to + 1 is the number of indices in the muhi-index iL + 2P — M . The 
parenthesis refer to the complete symmetrization of indices. 
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using A(l/ri) = —inSi (with ri = |x — yi|), Eqs. (14. 5 P then yield: 

= '^Si'vldj (^^^ +li^2 + 0(5) , (4.6a) 
Vf = ^S'^d, fl^ + 1 o 2 + 0{3) , (4.6b) 



2c ' ' \ri 

Wf^ = J± [sfvi^ - b'^S'tv'^ dk (^^^+1^2 + 0{3) . (4.6c) 

Notice that the leading-order of the spin part of these potentials has the typical dipolar 
structure d{l/r), while it is simply 1/r (monopolar) for the non-spin part. Note also the 
important fact that starts at 0{3) and has no 0{1) contribution. This makes the 
structure of the spin part of potentials different from the one of the non-spin part, and shifts 
a number of contributions to higher PN order. In particular, the non-compact support part 
of Wfj, coming from n^(—diVdjV), starts only at C(3) and not 0{1). This in turn implies 

that the leading order of Wfj has the simpler structure of a compact-supported term, which 
makes a number of sources of higher potentials less non-linear, such as the non-compact- 
supported term D:^{WijdfjV) in X^. 



C. Higher order compact-support terms 

We turn now to the higher-order PN computation of the compact-supported parts of 
potentials, i.e. whose sources include a, ai or cTjj. We will take the potentials as an 
example and give only a schematic view of these computations, since their implementation 
is relatively straightforward. By expliciting the expression (12. 8p of the stress-energy tensor, 
we get the following structure for the compact-supported sources: 

11 1 . 

ai(x,t) = ftiuSi + UibSi + dt{i^mSi) + di{ulj^Si) , (4.7) 

-9W v-gi^) v-^w 



indifferently for bodies 1 and 2, where the subscripts M and D refer to the monopolar and 
dipolar parts of the stress-energy tensor. The quantity fliu already intervenes in the study 
of the equations of motion without spin, and its 3.5PN expression in terms of the metric 
potentials is given in Eq. (4.3) of Ref. j58|. Next, we Taylor expand the retardations (14. ip 
for each term. The time derivatives are pulled out of the integrals, which are then evaluated 
using the rules (14. 2 p for the delta functions and their derivatives, keeping in mind that the 
factors 1 / y/—g in (14. 7p depend on the field point x. 

As a check of the method, we also computed all the compact-support parts of the spin 
potentials using a fully non-distributive prescription, i.e. keeping the dependence in the field 
point x as long as possible (also when evaluating the quantities such as i^m), and taking the 
Hadamard partie finie at the very end of the calculation. We obtained no difference with 
the pure Hadamard distributive prescription (13. 8p . 

Notice that , for instance, gets contributions from both the dipolar and monopolar 
parts of the stress-energy tensor. For the dipolar part, the quantities /iio, i^m and ul^) 
display an explicit spin factor. For the monopolar part, spin terms appear indirectly in two 
ways: through the spin contributions to the metric potentials in the expression of fliu, and 
through the acceleration replacements when evaluating time derivatives. 
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D. Non-compact support terms 



The non-compact supported potentials are more complicated, since their source is itself 
made of potentials and does not allow the simple integration procedure with delta functions. 
However, if the source has a simple structure, some integrations can be computed explicitly to 
obtain a solution valid for a generic field point x. One encounters often the source structure 
9(l/r)(9^(l/r), which arises at leading-order in all dA^^dB^ terms, where A and B are 
chosen among {V, V^j} and {V, ViWij} respectively, as can be checked from the leading-order 



structure fl4.6p . Using the same techniques as in previous works 
using the function g such that 



e.g. 



53 



in particular 



^9 



rir2 

ln(ri +r2 + ri2) 



we find the following relations: 



ri 



a: ( - I d,k ( - 



ri 



1 
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d,,k In(ri) + + d^'d, - 6^''d.:} (^^^ 



-dld%g 



(4.8a) 
(4.8b) 

(4.9a) 
(4.9b) 



where for instance d] = d/dy\. With these relations in hand, completing the results for the 
compact-support parts, we are able to straightforwardly compute the leading order 0{1) of 
, Rf, as well as the IPN relative order C(3) of Wfj, as indicated in Table H 
We now turn to the computation of the partie finie of Poisson integrals, of the type of 
the first and third term in fl4.ip FI When direct integration is not possible (at least using 
known results) for the non-compact support part of non-linear potentials, one by-passes 
the difficulty by directly computing the Hadamard regularized value of the integral at the 
points yi or y2. By doing so, one will not have access to the full information about the 
potentials (hence the metric) everywhere in space, but only the information relevant for 
the computation of the equations of motion. Since the computation is different for each 
derivative of the potential, one has to figure out which derivatives of which potential will be 
needed in the final computation. Table |T] already gave the required spin parts of potentials. 
Furthermore we find that two new non-spin (NS) potentials, featuring two derivatives, have 
to be computed in addition to known previous results (53| : 

(djkRf^)^ to order 0(2) , and (djkY/^^)^ to order C(0) . (4.10) 

Following Ref. jH^I (to which we refer for details and definitions), we define generic Poisson 
integrals and twice-iterated Poisson integrals: 



P(x) 
Q(x) 



.±Pf 
47r 



d^x' 



-Pf 



|x — X' 

dVlx 



x'|F(xO 



(4.11a) 
(4.11b) 



^ Pure spatial integrals like the second term of (|4.ip , where the integrand is independent of x, are computed 
using the Hadamard prescription for singular integrals, following the procedure described in Ref. 54 1. 
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where the source F is singular at yi and y2. The symbol Pf stands for the Hadamard partie 
finie for integrals, which depends on two arbitrary constants Si and So- The results for 
the extended Hadamard partie finie (in the sense of Eqs. (5.3)-(5.4) in jSJ]) of the Poisson 
potentials P and Q are: 



(^)l = 


-Ipf 

47r 


iQ)i = 


4:71 


(5.P)i = 


-ipt 

471 




4ir J 



d^^ 



ri 



-F(x) + 



Inf^ 

Sl 



[rlF) 



1 ) 



d=^x^F(x) + In 



ln(^ 



(rinlF)i , 



(4.12a) 
(4.12b) 
(4.12c) 
(4.12d) 



In these formulae, the first term would correspond to the naive prescription of taking directly 
x = yi inside the integrals (14. lip , and the second term features a regularization constant 
ln(r'^) coming from the singular limit x' — >■ yi. Notice that the constant Si automatically 
cancels between the two terms, while the Pf term may induce a dependence on so that 
the result depends on two constants and S2- For the present work we had to extend these 
formulae to the case of a double gradient, in order to compute the potentials f l4.10p . Using 
the same method as in Ref . 5J] , we obtained the corresponding required expressions for the 
double derivative: 



(dijQ)i 



471 



47r 



Pf / d'^x 



d'^x 



F(x) + ln( -i) ((3<-5^^)F)^ 



(F)i , (4.13a) 



■F(x) + 



In 



ij 



n 



(4.13b) 



One readily checks that these expressions yield the correct results when contracted with 6"^^: 
namely (AP)i = {F)i and {AQ)i = 2{P)i. Crucial for this check, notice the last term in 
Eq. fl4.13al) which stems from a distributional contribution obtained when evaluating the 
two derivatives under the integral fl4.11al) . It can be absorbed by replacing, in the first term 
of fl4.13al) . the ordinary factor (3n'*/ — 6^^)/r'^ by the distributional derivative dlj{l/r[). 

With these tools in hand, we were able to compute all the needed non-compact support 
parts of potentials. Importantly, we get no contribution of the singular parts proportional 
to ln(r'^) in Eqs. (14.120 . in any of the spin parts of potentials listed in Table HI This gives a 
strong indication that the Hadamard regularization (actually the pure Hadamard-Schwartz 
version of it) is sufficient to deal with all spin parts of potentials. This is not surprising 
because the calculation we are doing is only of 2PN relative order, while past experience in 
the non-spin case 53|, |55| says that Hadamard's regularization fails no earlier than at 3PN 
relative order. 



V. DIMENSIONAL REGULARIZATION FOR ONE POTENTIAL 

By contrast, we have found that in the Papapetrou part of the acceleration, which is 
specific to the spin-case (and outside past experience), we do get a contribution proportional 
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to \n{r[) in one of the new evaluations we had to perform, namely when computing the 
double-derivative potential {djdkY/^'^)i at Newtonian order using Eq. fl4.13ap . All the other 
NS potentials, including the other new one {djdkRf^)i we had to compute, could be safely 
obtained with the Hadamard regularization of Sec. IIVI 

The appearance of this In(r^) tells us, on the contrary, that Hadamard's regularization is 
insufficient for the potential {djdkY/^^)i, and that we a priori need dimensional regulariza- 
tion in order to get rid of possible ambiguities in the final equations of motion, as was the 



case of the non-spin 3PN equations of motion in |45l . |53|. On the other hand, when consid- 
ering directly the Papapetrou contribution to the acceleration, we see that the dangerous 
contribution exists in only one term, with index structure S('^{didjYj^^)i, and we find that 
the problematic logarithms In(r^) cancel out in the final result because of the antisymmetry 
of Si . Therefore, we expect beforehand the pole oc {d — 3)~^ we shall obtain in dimensional 
regularization to actually vanish in the final acceleration. Nevertheless, even though the 
result will be pole-free, we know that the finite part oc — 3)° should play a crucial role, 
and therefore we must a priori apply the powerful but tedious procedure of dimensional 
regularization to the problematic potential {djdkY/^^)i. 

We refer to jiH for precise definitions and technical details about the method, which 
consists in studying the problem in d spatial dimensions, treating d as a complex variable, 
and taking the analytical continuation of the obtained results when d — )■ 3. We will associate 
a superscript (d) to quantites which are defined in this d dimensional setting. One consid- 
ers a class of functions regular everywhere except at the points yi and y2, and admitting 
expansions of the type (VA^ G N): 

F('^)(x)= rr'-f^lM + oir^), (5.1) 



and similarly around y2, with e = d — 3, and Qq, qi, Pq G Z. As long as there are no terms 
in these expansions with both p < and q = 0, and that there is no angular dependence of 
the p = 0,q = term in the expansion, which is always the case in practice, the presence of 
the factor rf^ allows one to write, by application of the analytical continuation on e: 

(i.WGW)(yO = F('^)(yOG(''Hyi) , (5.2a) 

FW(x)5W(x-yi) = F('^)(yi)5^'nx-yi) , (5.2b) 

which justifies the use of distributivity in our previous Hadamard three-dimensional com- 



putations. Indeed, it was checked in previous work j45| on the 3PN non-spin equations of 
motion that the pure Hadamard prescription agreed with the d ^ 3 limit of the dimensional 
prescription for all the terms except those developing a pole. As explained in Sec. IIV C\ 
although we did not perform a full dimensional regularization computation of all quantities 
to be evaluated at yi or y2, we checked that distributive and fully non-distributive prescrip- 
tions for the computation of compact-supported potentials yielded the same final result, 
which gives us confidence in this distributive prescription]^ 

The treatment of Poisson integrals such as f l4.1ip goes as follows. Defining P^^^ (x) as: 



^ Within the dimensional regularization scheme, distributional derivatives can be treated as Schwartz dis- 
tributional derivatives, e.g. using Gel'Fand-Shilov formulae in d dimensions. 
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where k = T (^^) /vr^~, analytical continuation allows us to take derivatives under the 
integral, to set apart the distributional contribution, and to make directly the replacement 
X — )■ yi under the remaining integral. We obtain: 

(a,,pW)(yO = - 2) I d>^^ ^<-'^^V W(x) + ^i^^'^Hyi) , (5.4) 

where the second term on the right-hand side corresponds to the distributional contribution. 
It reads -F^'^''(yi) = iJq^q, according to the expansion (15. ip . and must not have any angu- 
lar dependence (since dimensional regularization does not include angular averaging in its 
definition, unlike Hadamard's partie finie). We set here to all the terms ff^iflfqini) with 
g 7^ since they are cancelled by analytical continuation (and have in general an angular 
dependence). We also assume that i/qq^ has no pole in e which is true, since no poles are 
generated at the level of the sources of potentials for this calculation. 

When d — )■ 3 we have F^^^ — )■ F, where F is the corresponding function in the Hadamard 
3-dimensional context. That function admits the expansion 

F(x)= Yl rlf,{n,) + o{r^). (5.5) 

The Hadamard partie finie of the function is defined by the angular average over the unit 
direction iii of the term p = 0, say 

(F)i = (/o(nO>. (5.6) 

The 3-dimensional coefficients i/p(ni) are related to the e — )■ limits of the d-dimensional 
coefficients i/pfg(ni) appearing in Eq. (15. ip by 

fp{ni)= Yl f?>^h (5.7) 
1 1 

qo^q^qi 

which in turn gives us the link between the analytical continuation when c? — )■ 3 of F^'^\yi) 
and the Hadamard partie finie {F)i: 

(F), = F'~'\y,) + J2 (/S(^i)>- (5-8) 

Next, we address the link between the Hadamard- regularized Poisson integral (14. lip and 
its dimensional regularization version (15. 3p . Defining: 

V{d,,P){l) ^ (9,,pW)(yi) - {dijP), , (5.9) 

and following the same steps as in Ref. j45|, we obtain, working at the zeroth order in e, 

V{d,,P){l) = Y (- + ^lnri)((rfni^-<5^^)/g> 

qo^q^qi ^ ^ 



17 



E + (5-10) 

The result of the PN iteration of the metric, starting with the (i + 1-dimensional Einstein 
equations, and the corresponding definitions for the metric potentials are given in Sec. II of 
Ref. [isj. In particular, the d- dimensional definition of the dangerous potential Yi is given 
by (2.12f) there]] We applied the previous method to the computation of the quantity 
'D[dj]^i){\) defined above. This required computing the singular expansion of the needed 
non-compact support sources in d dimensions up to the order p = 0, which was done following 



the methods of Ref. j45|. Our result for the pole part is quite compact. 



V{d,,Ym = \^-^v[,d]^u + 0{e') . (5.11) 

Since this expression is symmetric by exchange of the two indices i and j (or /c), it doesn't 
contribute to the final equations of motion where only the contraction S(^{dijYi^^)i is in- 
volved. This confirms that the final equations of motion in dimensional regularization are 
directly pole-free. Furthermore, we also checked that the finite part oc coming from 
Eq. (15.1 op also cancels out in the final result {i.e. after contraction with the spin tensor 
Si^). Although the latter fact could not a priori be guessed beforehand, it implies that 
finally the dimensional regularization was superfiuous for the computation of the dangerous 
term^f ((9ijF;f^)i. This result gives us further confidence that the pure Hadamard-Schwartz 
regularization is sufficient for all our calculations. 



VI. RESULTS FOR THE EQUATIONS OF MOTION AND PRECESSION 

Since higher-order results take the form of quite long formulae with many similar terms, 
we tried to adopt a systematic presentation. First, we split the expressions in terms of powers 
of G and according to the powers of the masses mi and m2. Each spin-dependent term is then 
denoted using the convention that gathers all the terms featuring G^^^m\m\. For any 
two vectors a and b, we use the notation {ah) for the scalar product, i.e. {ah) = a - b = a* 6*, 
and we define {Sah) = S'^^aV . 

Before presenting the results, let us address an important point: every time an expression 
has at least one free index, it might admit several equivalent writings, in terms of the vectors 
appearing in the problem and of the spin tensors. Indeed, any of the final results, which are 
functions of time only (such as the acceleration of one of the bodies or the time derivative 
of one of the spins), can be written with the three vectors n^g, ^i, and with at most one 
occurence of either one of the spin tensors and S2 since we are working at linear order 
in spins. As the spatial indices run on three different values, we have the two identities: 

S^'^a^^^ajhkCi = , (6.1a) 
amS"'^'a^h''c^^ajhkCi = , (6.1b) 

^ Notice that because we have already seen that the pole 1/e will necessarily cancel out in the final result, 
we do not have to worry about extending the Papapctrou part of the force p.lOp to d dimensions. I.e., 
we can evaluate the term djkYi assuming it carries the 3-dimensional coefficient, which is +8 in this case, 
see one of the terms in Eq. (|3.f Ocp . 
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where the brackets indicate antisymmetrization over the indices, where S^^ is one of the 
spin tensors and where (a, b, c) is a permutation of the set of three vectors {nu, Vi, ^2)- The 
number of identities of this kind that one must take into account depends on the number of 
vectors and tensors at disposal, and on the number of free indices. As an example, ( 16 .lap 
with S = Si, a = ni2, b = Vi and c = f 2 gives, once expanded: 



Hence, one must keep in mind that there is no unique writing for the results we are going 
to present, and take this into account when comparing two expressions. This becomes par- 
ticularly important when using the method of undetermined coefficients (for instance when 
looking for a contact transformation between harmonic and ADM variables, see Sec. IVII D|) : 
the system of independent equations that these coefficients have to solve is to be determined 
only after taking into account the complete list of these identities. Notice also that the 
use of an antisymmetric spin tensor reduces the number of these identities compared to the 
use of a vector spin variable and a Levi-Civita symbol, which was one of our motivations 
for changing the spin variable with respect to Paper I. A straightforward but cumbersome 
work around of this problem when comparing results is to project everything in an arbitrary 
orthonormal basis. 

A. Spin evolution equation 

First, we give the spin evolution equation, i.e. the equation giving the time derivative 
of the spatial components of the spin tensor up to 2PN order. This result was in fact 
already contained in Paper I, under the form of a precession equation for the spin vector 
SpBB; and it is a mere matter of traduction between the spin variables. Beware that our 
spin tensor is not of conserved norm, S^^S^^ 7^ const beyond leading order. This equation is 
needed each time we perform order-by-order reduction when evaluating time derivatives. 

Since we are working at 2PN relative order, it is sufficient to know the 2PN or 0(4) spin 
evolution equation. This also means that the amount of non-linearity in this computation 
is less that in the computation of the acceleration. Indeed, expliciting Eq. f l2.6ap in terms of 
potentials, we see that we need only 2PN potentials, and we can ignore the spin contributions 
in these potentials, working at linear order in spins. However, notice that, since the leading- 
order spin contributions to the total angular momentum of the system is of the form Si/c + 
S2/C, with Si^2 spin vectors, the order of this spin evolution equation required for finding a 
conserved total angular momentum at 3.5PN order is not 2PN but 3PN. 







- n\2vl{Sini2V2) + n\2{ni2Vi){SiViV2) - v\{Sini2Vi){ni2V2) 
+ vl{ni2Vi){Sini2V2) - v\{SiViV2) + S'^ v{{ni2Vi){ni2V2) 




(6.2) 



Defining 



^2B%^ + ^^2^N + O (6) 



(6.3) 



dt 



we obtain for the IPN order: 



Gm2 



2S'^ {m2Vi2) + ^n\2S'M: 



(6.4) 



IPN ~ 
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and for the 2PN order, splitting the result in terms of powers of G and occurence of the 
masses as explained above: 

G G^ 

^2PN = —Pq,i^2 + — [/3^imim2 + Pq^2^V\ ) (6-5) 

where 

+ vtS^^n\, [-?,{n,,v^f - 4(t;i2t'2)] + 4'5p< [4(ni2t'i2) + 2(^12^2)] 

+ S^^ [-?,{ni2Vi2){ni2V2f + 2{ni2V2){vi2V2)\ , (6.6a) 

PI, = 32{n,2Vu)42Sf'n'l2 - 14nf25?^f2 - 12i;[^5pn^2 

- 24'5?'nj2 + [2{ni2Vi2) + 2(^21^2)] , (6.6b) 

/3g2 = -4(ni2i;i2)nf25?'^t2 + 2v^uS'M2 - 2{nuvu)sr • (6.6c) 

Note that, with the series of potentials already computed for the 3PN equations of motion 
without spins 45|, |59|, we are able to control the precession equations up to the next 3PN 
order; this will be investigated in future work. 



B. Acceleration 



The spin contributions in the acceleration have the following structure, with our FN 
counting valid for maximally spinning-objects: 



-^ = Al, + -A 
at ( 



IPN + — 4 1.5PN + — 



5 



A 



2PN ~l~ ^ - 

SS 



2PN 



A 



AXt,^ + A^ 



13PN -r ;^^3PN 



A\ 



3.5PN + 4*3-5PN 



0{8 



2.5PN + 4* 



2.5PN 



(6.7) 



where the S subscript indicates the spin-orbit contributions, and the SS subscript indicates 
contributions that are quadratic in the spins and which we neglect in this workjl For the 
leading-order spin contributions, we get 



^1 4*1-5PN 



G 



12 



[m2 {SSl'n{2{ni2Vi2) + 6{Sini2Vi2)n\2 - SSl'vl2) 
+mi (6S'2^n{2(?^i2^^i2) + 6(S'2ni2t;i2)ni2 - 45*2^^^12)] 



(6.8) 



V\ — Vn. 



It depends on the two velocities vl and only through the relative velocity v{2 
This is imposed by the Lorentz invariance of the harmonic-coordinate equations of motion, 



Notice that, when doing the calculation in the original harmonic coordinates, one actually obtains some 
3PN smn-orbit contributions. However, these are pure gauge, as was already explained in Appendix A of 



Ref. 



60[, and are eliminated by the gauge transformation 



SX'' with dX° = and 



6X' = - 



G2 



771182 — rn2Sl'' 



whose effect on the accelerations is Sa\ = 5a\ = d^SX^ /c\t^ + 0(8). This gauge transformation obviously 
respects the harmonicity condition in a perturbative sense, since \I\5X^ — 0(8). 
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which reduces to a Gahlean invariance at leading order (see Sec lVIIBj) . For the IPN relative 
order contributions, we get the following: 



G 

mi 4*2.5PN = — [«0,l"^2 + «l,o"^l] + — ["0,2^^2 + «l,l"^l"^2 + a2,o"^l] 



(6.9) 



12 



12 



where the coefficients are completely equivalent to those obtained in Paper I and reac^l 



"0 1 — '-'i 



'1 "12 



15 3 

Y (^12^12) ('^12^^2)^ - 2(^12^12)^12 - 3(r2l2t'l2)(^^12^^2) 



-3(^12^2) (^'I2f 2) - 7t('^12^'12)^'2 



9 3 3 

-3(ni2fl2)(?^12^^2) + -(?^12^2)^ + -fl2 + 6(t;i2^^2) + 



(6.10a) 



+ ('S'ini2'yi2)'^i2 [-15(rii2f2)^ - 3t'i2 - 12(t'i2t;2) - 3^2] 
+ iSini2Vi2)v{2 [-3(ni2t;i2) - 6(ni2t;2)] + 3(ni2t;i2)(5'i?ii2t;i2)^'2 
+ 3(ni2i;i2)(5'ini2f 2)^12 + 3(ni2i;i2)(5'ini2?;2)^^2 
- 3(5'it;i2t;2)^^l2 - 3(S'it;i2t;2)t'2 , 

"1,0 = 'S'2^'^12 [-15(ni2fl2)(?^12f2)^ - 6(rai2fl2)(t'l2f2) + 6(^12^2) (fl2^'2) " 3(r2i2f 12)^2] 

+ 5'2^^i2 [6(ni2t;2)^ + 4:{vi2V2) + 2V2] + iS2ni2Vi2)n\2 [-15(ni2f2)^ - 6(f 121^2) - 3^2] 
+ (52ni2fi2)fl2 [-6(7^12^12) - 6(rai2f2)] + 6(ni2t;i2) (52^12^2)^12 

+ 6(^12^12) (5*2/112^2)^2 ~ 4(52t;i2^^2)^42 ~ 4(52^12^^2)^2 ' 



a, 



0,2 



a 



1,1 



= -6(r;,i2t'i2)5i-'?7,{2 + 65i-'?;i2 - 12(5ir2i2i'i2)'T-i2 

= -14(ni2t;i2)5i-'n-{2 + 145i-'t;i2 - 26(5ini2t;i2)?^ 

- 16(r?,i2t;i2)52-'r2{2 + '^'282 VI2 ~ 20(52^12^12)^12 5 
31 



12 



^2 = Sl''n-^ 



'2 "-12 



-(nuvu) + 2(ni2t;2) 



+ y52^i;^2 - y(52ni2Wi2)ril2 • 



(6.10b) 
(6.10c) 

(6.10d) 
(6.10e) 



Finally, the main result of our work, namely the next-to-next-to-leading 3.5PN spin-orbit 
contributions to the acceleration, reads: 

G G^ 

^1 43.5PN =-3- [7o,l"^2 + 7l,0"^l] + -r [7o,2"^2 + 7m"^1"^2 + 72,0^^?] 



12 



12 



G' 



where 



7n 1 = S'^n- 



12 



105 



+ — [7o,3"^2 + 7l,2"^l"^2 + 72,l"^?"^2 + 73,0"^?] 
^^12 



15 3 

(?^12t'l2)(ni2f2)'' + — (ni2t;i2)(?^12^^2)^t'i2 - -(ni2t'l2)f 



(6.11) 



12 



15 



ni2fl2)(?^12t'2)^(fl2t'2) 



15 



?^12f2)^(^'l2W2) - 7T('^12^12)fi2(t'l2t'2) 



The link between the spin variable used in Paper I and the present spin tensor is provided in Appendix Rl 
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-^('^12^^2)f?2(^12^^2) - ^ (^12^^12) K2^'2)^ " 3(ni2t;2) K2^^2)^ 



-^{ni2Vi2){ni2V2fvl - ^{ni2Vi2)vl2V2 



+ 7::{ni2V2){Vi2V2)vl - -(ni2Wl2)t'2 



•^1 ^12 



15 



45 



:(?^12fl2)(t'l2f2)f2 



9 



2„.2 
12 



Y('^12^^12)('^12^^2) -y('^12^^2) + 2 (^12^12) ('^12^^2)^^12 " ^ ('^12^^2) 

3 3 

+ g^l2 + 3(ni2t;i2)(^12^^2)K2^'2) " 12 (ni2f 2)^ K2^^2) " 2(^12^2)^ 
-^(?^12^12)(^12^^2)^2 + ^('^12^'2)^f2 + ^^12^2 + 3(t;i2^2)^'2 + ^^2 



{Sini2Vl2)n\2 



105 15 3 

— (^12^2)'^ + Y(?^12f2)^Wi2 - -^^U + 30(r2i2f2)^(^^12^^2) + 3(t;i2f2)' 



15 

y 



-^{ni2V2fvl - ^t'i2f2 - Q{Vl2V2)vl - 7W2 



3 

-1 

4 



+ ('S'ini2t;i2)i' 



12 



+ (5'ini2fl2)t^2 



(S'ini2f2)t^ 



12 



15 3 

Y(?^12^12)(^12f2)^ + 15(ni2f2)^ + -(?^12^'12)^'?2 + 3(ni2f2)t'?2 

9 

+3(rai2fl2)(fl2f2) + 9(7212^2) (^^12^^2) - -('^12^12)^2 - 3(ni2f2)^'2 

15 3 

-Y('^12fl2)('^12^^2)^ - 2('^12^12)^^12 " 3(rai2t'l2) (f 12^^2) 



-3(ni2f2)(^^12f2) - 7t('^12^12)^^2 



15 3 

— (?^i2^^i2)(?^i2^^2)^ - -j{ni2Vi2)vl2 - 3(ni2t;i2) (^^l2^'2) 



+ {Sini2V2)vl 



+?>{ni2V2){vi2V2) - - {ni2Vi2)vl 



15 3 

-y ('^12^^12)('^12^^2)^ - -('^12^^12)^'l2 " 3(^12^12) K2^^2) 



-3(^12^2) (t'l2f 2) - -('^12^^12)^^2 



9 3 3 

+ {SlVi2V2)v\2 -3(r2i2t'l2)(ni2t'2) + -(?^12t'2)^ + -f?2 + 6(fl2t'2) + 

9 3 3 

-3(ni2fl2)(ni2f2) + -('^12^^2)^ + -VI2 + '0{Vi2V2) + 



{SiVi2V2)v2 

105 



(6.12a) 



11,0 — ^2 "'12 



{ni2Vi2){ni2V2)'^ + 15(ni2t;i2)(ni2t'2)^(^^i2W2) - 15(^12^2)^(^12^2) 
15 

6(ni2t;2)(^^12^^2)^ - ^ (^12^12) ('^12^^2)^^^2 - 3(ni2t;i2)(^^12f2)^^2 



+3(ni2^;2)(^^12^^2)f2 - ^('^12^^12)^'2 



22 



12 



15 



{ni2V2) - 6(ni2t'2) {V12V2) + 3(ni2'y2) + 2('Ui2f2)'y2 + o"": 



+ (5'2ni2fl2)?^ 



12 



105 



(ni2t;2)'^ + 15(ni2t;2)^(^i2^2) 



15 
Y 



{ni2V2fvl 



3(t'l2^2)^2 ^ 7^'2 



+ 
+ 
+ 



70,2 



+ (5'2ni2fi2)t^l2 [I5(ni2^;i2)(ni2^^2)^ + 15(ni2t;2)^ - 3(ni2t;i2)^^2 " 3(ni2t;2)^^2] 

+ (5'2ni2t;2)^^l2 [-15(ni2t;i2)(?^12^^2)^ - 6(ni2t;i2)K2^^2) + 6(ni2t'2)K2^'2) - 'i{'^l2Vl2)vl] 

{S2ni2V2)v\ [-15(ni2fi2)(?^i2t'2)^ - 6(rai2?;i2)(t'i2^'2) + 6(ni2t;2)(fi2^^2) - 3(^12^12)^2] 
(S'2t;i2f 2)^12 [6(^12^2)^ + 4(t;i2t'2) + 2^2] 

(^2Wl2^^2)^^2 Hni2V2f + 4.{V12V2) + 2^72^] , (6.12b) 

S'^n{^ [I8(r2i2fi2)(^i2^^2)^ - 3(^12^12)^12 + 6(ni2t;i2)(t'i2f2) 

-6(r2i2f2)(^'l2f2) + 3(^12^12)^2] 

Siv{2 [-6(ni2t'i2)(rii2f2) - 12(ni2W2)^ + 3^12 - 3^2] 

(S'ini2fi2)ril2 [36(rai2f2)^ - 61^ + 61^2] + ('S'ini2fi2)fl2 [-6(rii2t'i2) + 12(ni2W2)] 
6(ni2fi2)(5'ini2t^i2)^^2 - 6(^12^712) (5'i?ii2t;2)^42 - 6(ni2t;i2)(5'ini2t;2)^^2 

6(5'i'Ui2f2)'yl2 + 6(S'it'i2f2)f2 , (6.12c) 



+ 
+ 
+ 
+ 



7i, 



Sin{2 



375 



(7212^12)^ - 33(ni2t;i2)^(ni2t^2) + 45(ni2t;i2)(^i2^^2)^ 



177 



(7212^12)^^12 



-14(rai2fl2)(t'l2f2) - 14(ni2t;2)(t'i2f2) + 7(rai2f 12)^2] 



12 



271 



+ 

+ (5'irai2t'l2)n^ 



(r2i2t'l2)^ + 39(72121^12) (?^12^'2) - 32(7Zi2t^2)' 



73 



'12 



12 



663 



'^12^'12)^ - 42(7112^12) (?T-12^'2) + 90(r7i2f2)^ + 



28(1^12^2) - 7u, 
135 



+ 



+52(vi2t;2) + 131-2] 

(S'i7ii2fl2)t^l2 [122(7212^712) + 20(7ii2t;2)] - 14(71121^12) (5'l72i2i;i2)^72 

14(77,12112) (5'l71i2f2)f 12 - 14(r7i2ll2) (5*1721212)12 + 14(S'ili2l'2)ll2 + 14(5'ili2l2)l'2 

n\^2 [43(7212112) (''^1212)^ + 16(7212112) (112I2) ~ 16(7^1212) (112I2) + 8(72,12112)12] 

5*^1^2 [-24(721212)' - 12(11212) - 6I2'] 

(^2ni2ii2X [60(721212)' + 20(11212) + IO12'] 

20(2212172) (5'272i2i;i2)^^l2 - 16(72i2l'l2)(5'272i2l'2)^42 " 16(r2i2f 12) ('S'272 i2272)^'2 
12(5'2ll2f2)fl2 + 12(S'2fl2f2)t'2 , 



(6.12d) 



72,0 



S2'^12 



1815 



801 



(7212112)^ - 51(7212112)^(^1212) + 54(7212112) (^1212)' - 6(721212)' 



(7212112)1 



12 



31 

12(721212)112 + y(?^12112)(11212) 



39 



(72i2l2)(ll2l2) 



31 

+ -r('^12ll2)l2 - ('^12l'2)l2 



4 



12 



1087 



('T'12^'12)^ 



59 



56(7212112) (rii2i2) - y (^^■12^2)' 



269 



''12 



23 



+ {S2ni2Vi2)n\2 



1797 



93 



323 



8 



12 



+ {S2ni2Vl2)v 



12 



+ {S2ni2V2)vl2 



,45 
+ Tr(^i2^2) + -r^: 



67(^12^12) + Y(rii2f2) 
31 

— ^(^^12^^12) + '2{ni2V2) 



24(ni2f 12) (5'2ni2f 12)^2 
31 



(S'2r^l2^;2)^'2 



(^^12^^12) + '2{ni2V2) 



23 23 

+ Y(5'2t'l2^'2)^'l2 + Y(5'2fl2W2)w2 ) 



76,3 



15 

227 



'^12^^12)5'i^'^' 



12 



15 



Sl^vi^ + 15(S'ir;.i2yi2)n' 



12 ' 



309 



7l,2 = ^ (^12^^12) ^l^«i2 - ^-^l^^' 



691 



12 



{Sini2Vi2)n 



12 



32{ni2Vu)S'2'n{^ - 24^2^ 17^2 + 42(52ni2t^i2)ni2 

79 • 63 

y(ni2t;i2)^/n{2 - y 



72,1 = — ('^i2^^i2)5'l^'T'i2 - —Si^l^ + 58(S'ini2t;i2)nl2 

615^^t;j2 + 



2 "-12 



251 



(jT-12^'12) - 14(ni2t;2) 



257 



{S2ni2Vi2)n 



12 ) 



73,0 ~ ^2 '^12 



73 



rai2fl2) - 14(?2i2t'2) 



119 



343 



'S'2''^'i2 + -ir{S2ni2Vi2)n 



12 



(6.12e) 
(6.12f) 

(6.12g) 

(6.12h) 
(6.12i) 



VII. CHECKS OF THE RESULTS 



A. Conserved Energy 



An important feature of all the spin-orbit contributions we have computed, is that they 
are associated with the conservative part of the dynamics, i.e. obtained when neglecting the 
radiation-reaction dissipative terms associated with gravitational radiation. Therefore these 
contributions should allow for the existence of a set of conserved quantities, namely energy, 
angular momentum, linear momentum, and center-of-mass integral. Here we have checked 
that the new terms computed in the equations of motion in Sec. IVI Bl admit corresponding 
contributions in the conserved energy, which take the following structure: 



E — Ej<{ + — EipN + — Ei.sPN + 



+ 



1 



s 



-E'sPN + E 3PN 

ss 



-E'2PN + -E'2PN 
SS 



H — 7 E2.5PN 



s 



+ ^ E3.5PN + o{- 

c' s 



The leading order spin-orbit contribution to the energy reads 

'^'"'(5ini2fi) + l 02 



-E'l.SPN 

s 



(7.1) 



(7.2) 



24 



while the 2.5PN result (IPN relative) was given in Paper I, and reads, once translated into 
the spin tensor variable S''-': 



-E'2.5PN 

s 



r 



G G'^ 

— eo,im2 + — [eo,2"^2 + ei,i''^i"^2] +1^2 



12 



r 



(7.3) 



12 



where 



eo,i = iSini2Vi2) 
+ (5'ini2f2) 



9 

3(7^12^^12) ('^12^^2) + 2(^12^2)^ + (^^12^^2) 



-3(ni2i;i2)^ - 6(ni2t;i2)(rii2^^2) - -(?^i2^^2)^ + v 



12 



+ {S1V12V2) [(r2l2t'l2) + 3(ni2t'2)] 

eo,2 = 2(S'ini2t;i2) - {Sini2V2) , 
ei,i = -2(S'ini2t;i2) - 2{Sini2V2) ■ 



(7.4a) 
(7.4b) 
(7.4c) 



The result for the 3.5PN contribution to the conserved energy has been obtained by the 
method of unknown coefficients and successfully determines the following unique result: 



-E3.5PN = -^V0,l1^2 + ^ ['70,2"^2 + »7l,l"^l"^2] 



12 



12 



G 



— [^,3^2 + r]i,2'miml + r]2,imlm2] +1^2 

^^12 



(7.5) 



where 

'70,1 = {Sini2Vi2) 



iSini2V2) 



+ (5'it;i2f2) 



15 45 

^(ni2fl2)^(?^12t'2) - 15(ni2t'l2)^(ni2W2)^ - y (?^12t'l2) (?^12f 2)^ 

"-^-('^12^2)^ + 7('^12^^12)('^12^^2)f?2 + 3(ni2t;2) V2 - 7('^12^^12)^(^^12f2) 

o 4 4 

9 1 

+3(ni2t'l2)(^12^2)(fl2^'2) + -('^12^^2)^(fl2^'2) + -^12(^12^'2) + (fl2^^2)^ 
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+3(ni2fl2)^^^2 + 12(ni2fl2)(ni2^;2)^^2 + y('^12^'2)^^^2 ^ ^12^2 + 2(t;i2^'2)^^i 

"15 45 

^(^^12^^12)^ + ^5{ni2Vi2f{ni2V2) + y (ni2t;i2)^(rii2^'2)^ 

15 9 

+ 15(ni2Wl2)(ni2t^2)^ + -^(?^12^^2)^ - 7t('^12^^12)^^^12 " 9(ni2t;i2) (^^12^^2)^'l2 

o / 

9 3 

--{ni2V2fvl^ + -17^2 - 6(ni2t;i2)^(^^12f2) - 12(ni2t'l2) (?^12^^2) (^^12^^2) 

-6(rii2^2)^(^^i2^2) + 2i;i2(^i2^'2) + (^^12^2)^ - 3(ni2t;i2)^f2 

-6(7212^12) (ni2t'2)^2 - ^{.nx2V2fvl + VI2VI 

3 9 9 

-^(^^12^^12)^ - 3(ni2t;i2)^(ni2t;2) - 2('^12^12)('^12^^2)^ - 2('^12^2)^ 

3 

+ ^(ni2t'l2)^^?2 + ('^12^^2)^^?2 + ('^12^^12)(^^12^^2) + (?^12^'2) (^^12^^2) 
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V0,2 = {SlUuVu 

+ {Sini2V2) 

+ {S1V12V2) 

Via = {Sini2Vi2 
+ iSini2V2) 
+ (5'it;i2f2) 



+ (^12^12)^2 + 3(rai2f2)^^2] 
53 



(7.6a) 



23(^12^12) (rii2f 2) 



109 



53 



35 



V12V2) - V2 
65 



('^12^'l2) + —('^12^^12) («12^^2) + 3(ni2W2) + —V 



12 



^7 29 

— ('^12^'l2) - y(ni2^^2) 



(7.6b) 



105 



141 



7 



('^12^12 



(^12^12) ('^12^'2) - 22(^12^2) + -^'12 + 18(^'12^2) 



161 



4 
123 



13 

^^12^^12)^ + 22(ni2fi2)(?^i2^^2) + 4(rai2f2)^ - y^L 



(^12^12) - 14(ni2t;2) 



?7o,3 = 2{Sini2Vi2) - -{Sini2V2) , 
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^1,2 = y (5'ini2t;i2) + -^iSini2V2) 



(7.6c) 

(7.6d) 
(7.6e) 
(7.6f) 



^72,1 = y (5'ini2t'i2) + ^{Sini2V2) ■ 

We leave for future work the study of the other conserved quantities at the same order, 
namely the total angular momentum, linear momentum, and center-of-mass integral. 

B. Lorentz invar iance 

Since we are working with the harmonic gauge condition which is manifestly Lorentz 
invariant, the global Lorentz invariance must be preserved by our calculations and must be 
manifest on our final equations of motion^ To check this, we follow mostly the presentation 
of Ref. 6l| and of the Appendix A of Paper I, with the difference that we are using a different 
spin variable, which will simplify this calculation. 

Let us consider two different frames (J-") and {J^'), the latter being related to the former 
by a boost of velocity V. The coordinates of a given space-time event P are in the 
original frame (J-") and x'^ in the boosted frame {J^'), both being related by the Lorentz 
transformation x'^ = A'^^.x'^, with: 

(7.7a) 
(7.7b) 

7 + 1 

with 7 = (1 — y^/c^)~^/^ the Lorentz factor. We denote the trajectories of the two bodies 
in the frame (J") by = {ct, yi) and y!^ = {ct, y2), and by y'l" = {ct', y[) and = {cf, y2) 



A° 




= 7 , 





= A°.= 


A*. 


= 6) + 



-7— , 

c 

72 v% 



The global Lorentz invariance is the one associated with the asymptotically Minkowskian space-time far 
away from the compact-support matter distribution. 
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in the boosted frame {J^')- The point is that we cannot compare them directly, because the 
simultaneity surfaces are different in the two different frames. To give a sense to simultaneity 
between the two frames, it is convenient to define an auxiliary event f2(ct, x), which can be 
choosen at will — it could be for instance the point of coordinates (ct,x = 0) in (J^), 
and whose coordinates in (J-"') will be (ct',x'). We define the two events Pi(ct, yi(t)) and 
P2(ct, y2(t)) on the two worldines, simultaneous to Vt in the frame (J-"), and similarly two 
events Qi{ct' ,y[{t')) and Q2{ct' .y^it')), simultaneous to VL in the frame {F')- ^ plays the 
role of an observer, for which the equations of motion evaluated on the simultaneity surfaces 
(f2,Pi,P2) and {^■,QiiQ2) must take the same fonctional form, as explained below. Next, 
we define the times ri and r2, such that the coordinates of Qi and Q2 in (J-") are (cri, yi(ri)) 
and (cr2, y2(T2)). We have, by construction, y'i^{i[) = A^^?/J'(ri) and similarly for 2. The 



link between y'i(t') and yi{t) is provided in Ref. [6l[. One obtains successively 




7 



7 + 1 



V 



(7.8a) 
(7.8b) 
(7.8c) 



where the right-hand sides are evaluated at t, and (Vri) means V- (x — yi). The expressions 
for the velocity and acceleration are obtained by taking the time derivative according to 
d'f. = 7(9t + '-fV^di. Similar expressions hold for 2, and transformations for quantities such as 
ri2 or are deduced in a perturbative sense from the first of these formulae. 

We turn now to the transformation rules for the spin tensor. As explained in Paper I, we 
have for any function f{t): 



fin) = fit) + J2 



n=l 



in— 1 



(7.9) 



which can be applied to the components Si"^ and Si of the spin tensor. Now, since these 
are the components of a contravariant tensor, we have 



S['\t') = A\Ksrin) , 
sfit') = A^^A\srin). 

Combining these expressions, we arrive at the desired transformation rule for S^: 



(7.10a) 
(7.10b) 



S'l'it') = SI 

+00 



1 C2 ' ^ 



n=l 



7 
c n\ 



in— 1 




27^ 1/M5f^,^27d^^, 



7 + 1 c2 dt 



dt 



(7.11) 



where S^^ can be further eliminated using the spin supplementary condition (12. 9p . This 
transformation is somewhat simpler than (A13)-(A15) in Appendix A of Paper I because 
the transformation for the spin vector S^-q^ is not as simple as in Eqs. ( I7.10p . 
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Lorentz invariance means that the acceleration in the boosted frame, obtained by trans- 
forming all variables according to Eqs. (17. 8p and (17. lip , when systematically truncated at 
the requested PN order, must take the same functional form as in the original non-boosted 
frame. We have verified that our result for the acceleration (16. lip passes this test. Note 
however that the test leaves a lot of freedom in the acceleration. In particular, one might 
add at the highest level 1 /c^ in the acceleration any quantity depending only on the relative 
velocity w ^2 = vl — V2 ^^^^^ P^^^ this test. 

Finally let us remark that the pure Hadamard Schwartz regularization which has been 
used for all the terms in this computation, finally yields equations of motion which are 
manifestly Lorentz invariant. This is another indication that there is no need for using a 
more sophisticated regularization such as dimensional regularization in this problem. 



C. Test-mass limit 

Another important check of our result is to take the test mass limit and show that we 
recover the equations of motion of a test particle orbiting in a black hole background. More 
precisely, we can check that: 

1. In the limit where one of the two bodies (say body 1) is a test particle with rrii — )■ 
and is spinless {Si = 0), its acceleration (16. 7p reduces to that of a spinless test particle 
orbiting around a Kerr black hole; 

2. In the limit where body 1 is a test particle with spin [i.e. mi — with constant ratio 
Si /mi), Eq. (16. 7p reduces to the equations of motion of a massless spinning particle 
orbiting around a Schwarzschild black hole. 

Note that in case 2. we need to work with a Schwarzschild black hole because we restricted 
ourselves to spin-orbit effects. In this section, it will be more convenient to rewrite the 
equation of motion (l3.5p -f lXB]) so as to make apparent the coordinate acceleration: 



dw^ d'^x^' „ „ ^v^'„r, \ 1 



v^'V ( -T\^ - ]+-A^'-^-^'] ^ (7-12) 



where we recall that J-"^ = —^^R^^^ S^"" and = 1/ ^y—gpo-vPv'^ /c 



1. Spinless test particle: equivalence with Kerr geodesies 

In the limit where mi — )■ and Si — ?■ 0, then = and S2 = const is a trivial 
solution of the equations of motion and the (spin part of the) acceleration of body 1 given 
by Eq. (16.111) . simply reduces to 

G G 
= -5-3- (6(ni2fi)v42 - ABl + &{S2ni2Vi)n\^ - Q^-^{ni2Vi){S2ni2Vi)v\ 
c ri2 c ri2 

q2 

+^T- (-16m2(ni2t'i)A*2 - 20m2(^2rii2t'i)<2 + 12m2S^) 
c '^12 

{22ml{ni2Vi)^2 + ^2ml{S2ni2Viyi2 - 2AmlBl) , (7.13) 



c'r 



12 
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where we have defined = ^[282 and B2 = v{S2 ■ For simphcity, we choose the origin 
of our coordinate system at the location of the central body 2 and suppose without loss of 
generality that the spin of body 2 points along the z axis: defining 5* = e^^^S'!^ , we impose 

Given the symmetry of the problem, it is of course more convenient to work with the 
spherical coordinates associated with our harmonic coordinates and with the associated co- 
ordinate basis that we denote {dr,dg,d^). In practice, we will show that our PN result 
recovers the dynamics of a test-mass in a Kerr background by comparing the explicit ex- 
pressions for the spin parts of the quantities f, 9 and (p obtained using Eq. fl7.13p and the 
geodesic equation in Kerr. In terms of r, 6 and (p, we have the following set of scalar quan- 
tities ri2 = r, (ni2t'i) = r, (>S'2ni2t'i) = — m2ar</) sin^ 6*, and vector components v{ = r, 

= 6, vf = 0, = 1; ^2 = — m2a/r, B2 = — m2ar0sin^^, B2 = —m2a<j) cos 6 sin 9 
and i?2 = ''Ti2a(^cot^ + r/r), all the other ones being zero. The components of the 
Cartesian acceleration are = f — rO"^ — r^^sin^^, = 6 + IrQjr — 0^ sin 6' cos 6' and 
a<^ = + 2r0 cot 6'/r + IQip cot 6. This yields 

, , Gm2a ■ 9 ^ Gm2a -9 9 ^ G'^mla ■ 9 ^ G^mla ■ 9 ^ ,^ , . 
(r)s = 2—^(f) sin^ 6 - 6^^0 gj^s q _ 8 ^ ^ (p sin^ 6 + 18 ^ ^ (p sin^ 6 , 7.14a 

Gm2a ■ ^ ^ Gm2a ■ ■ 9 ^ G^mia ■ 
(0)s = cos ^ sin ^ - 6^^f ^ sin^ ^ + 12 , cos ^ sin 6 

-2A^^4^(Pcos9sm9 , (7.14b) 



Gm2a , ■ ^ , Gm2a -n 9^ G'^mla 
—^UrO cot ^ - 2r - G^^r^^ sin^ ^ + —-^( 



G^m\a ■ 



+ 24 ' ^cotg . (7.14c) 

We would like to compare this with the equations of motion of a test particle in the back- 
ground of a Kerr black hole of mass vn2 and spin parameter a. Since our result is written 
in the harmonic gauge, we also need to work with the Kerr metric in harmonic coordinates 
rather than in the usual Boyer-Lindquist (BL) ones. A particular set of spatial harmonic 
coordinates (x^, and x^) constructed from the BL grid was obtained in Ref. j62[ and reads 



X + ix = (rBL - "^2 + ia) sin 6'bl exp i 

„3 



^>BL H In 



(7.15a) 



X = (rBL - ^2) cos 9bl , (7.15b) 



with r-i- = (m2 ± y m^ — a^). For the time coordinate, we can simply choose t = tBL- Here 
again, we will use the spherical coordinates associated to x^, x"^ and x^. The line element, 
expanded to linear order in the spin a, reads 

^ _r-rn2^^2 ^ '^_±^^^2 ^ ^ ^^)2 ^^^2 ^ ^^^2 g ^^2) 

r + m2 r — m2 



Note that in order to avoid heavy notations, we have set G = c = 1 in these non PN-expanded 
equations ( I7.15p - (l7.16p . Using ( I7.12p . in which we set to zero, for the coordinates 
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x'^ = {t, r, 6, 0), and developing at linear order in a and at 3.5PN order, we obtain expressions 
for r, 6 and (p that reduce to the results of Eqs. fl7.14p . but with the addition of the following 
extra pure gauge terms at 3PN order: 

{6r)g = -2—^-^r(j)sm^e , (7.17a) 
Se) = -2 „ f r(j) cose sine , (7.17b) 



s c"r'^ 



fs;p\ =_9!^^2r'^2rrecote-r\e' + <p'sm^e))-^^^ . (7.17c) 

The presence of these terms is simply a consequence of the fact that the coordinates fl7.15p 
are not exactly the same as the ones we used in our PN calculations: the harmonicity 
condition does not completely fix the gauge and we still have the freedom to perform a 
coordinate change — > + as long as dp{^J—ggP'^dai'^) = without violating the 
harmonicity condition. More precisely, the presence of these terms can be traced back to 
the metric element (7^0 which can easily be put to zero (modulo higher order PN corrections) 
by the coordinate shift 

(7-18) 

In terms of our harmonic coordinates, this translates into a shift = — ^"^^^^" ^ and 
^2 _ "^G^m^a ^ which can easily be seen to satisfy dp{y^—ggP^da-^^) = modulo higher than 
3.5PN corrections: since C,^ = 0{l/c^) (and it does not depend on time), this equation 
simply becomes A,^* = 0. 



3. 



2. Spinning test particle: equivalence with Papapetrou motion in Schwarzschild 



We now consider the limit mi — > while /mi is kept constant and we set S2 = 
since we only want to study linear effects in the spins (the spin-orbit terms involving 5*2 are 
precisely the ones that have been studied in the previous subsection). Here again, f 2 = is 
a solution of the equations of motion and the acceleration of body 1 reduces to 



(oi)^ 



G m2 



c^ri2 mi 
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m2 




mi 


G^ 


ml 


&r\2 


mi 
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m2 


CV?2 


mi 




ml 




mi 


G^ 


m\ 





"12 



-{ni2Vi)vlA\ - 'ivl{Sini2Vi)n\2 + -^vlB\ - 'i{ni2Vi){Sini2Vi)v\ 
[-&{ni2Vi)A'i - l2{Sini2Vi)n[2 + &B\) 



3 3 3 3 

--{ni2Vi)vlA\ - -vl{Sini2Vi)n\2 + ^^'i^i + 7:ii^i2Vi){Sini2Vi)vl 
o 4 o Z 



-3{ni2Vi)viA\ - 6vi{Sini2Vi)n\2 + SvfBl - Q{ni2Vi){Sini2Vi)vl) 



15 15 
y (7212^1)^1 + 15{Sini2Vi)n\2 - —B\ 



(7.19) 
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with A\ = ni2Si and Bl = vlSl^ . Choosing once again the origin of the coordinate system 
at the position of body 2 and moving to spherical coordinates, we obtain = 0, A{ = — S'[^, 
At = B{ = r^9S[' + r2_0sin2 9Sl'^, = -rS[' + r^sin^ OS^^, Bf = -f^[^ - r^esf'^ 

and (5'ini2fi) = r'^6S{^ + r'^^'^in^ 6S['^ . Plugging this into Eq. f l7.19p leads to the explicit 
expressions 



[r- s 



+ 



, G m2 
(?r m\ 
G m2 



, ml 



3 G m2 
2 (?r nil 



3 



6 



4M 



(^3f 



4 8 4 



2 -2 i2 • 2/1 

r r (h sm 



_G^m| 
c^r^ mi 



4/d2 <2 



15 ml 
2 c'r^ mi 



re 



= -3^^0sin2 esl^ - 3^^ ( rO-'Sl' + r6 
c^r mi &r mi 



sin 6 — —r 
> sin^ 0S'[ 
" ' sin^ 9Si 



^-^Sin^e 
(7.20a) 



—V 

2 



+ 6 



ml 
c^r^ mi 
mi 
c'r"^ mi 



3^^^5f 
c^r mi 



sin2<^ + 



3 G m2 



2 c''r mi 



1 



.sin' ^^['^ + ^^'0 sin' 



ism 



15 G^ m 



2 c'r'^ mi 



)sm 



(7.20b) 



6 



-6 



G^ mj 
c^r^ mi 
_G^m| 
cV^ mi 



^ G m2 
c^r mi 
3 G m2 



r^05[' + sin' ^^['^ + -Ov'Sl' 



1 . 



+ -^—v' r^0^r + ^0 sin' ^^i 
2 c'r mi \ 



4 ^ 



15 G3 



mr 



2 c'r^ mi 



(7.20c) 



where we have used the shorthand notation v"^ = (r^ + r^^^ + r^0^ sin^ 0). 

We want now to compare these results to their counterparts as given by the Papapetrou 
equation (]2.6bp written in the Schwarzschild background and in harmonic coordinates. A 
set of such coordinates and the corresponding form of the metric can be readily obtained by 
setting a = in Eqs. fl7.15p and (17.161) . We use (17.121) . including this time the Papapetrou 
force J-''^, and we obtain expressions for f, 9 and which, when expanded at the 3.5PN 
order, take exactly the form of Eqs. (17.201) . 



D. Equivalence to ADM results 



In this section, we compare our results to the ones obtained in Ref. [38| by a completely 
different method, using a reduced Hamiltonian formalism in ADM-type coordinates. In the 
following, we will denote all ADM variables with an overline. We follow the method of (32j . 
which compared to the harmonic-coordinate results in Paper I at the next-to-leading order 
C(5). In this formalism, the canonical structure is defined for the variables (positions 
of the two bodies), (canonical momenta) and Sa (canonical spins). Notice that these 
canonical spins are of conserved norm. The Poisson brackets for these variables read: 



{K.f,]=5a,5''. (7.21a) 
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\s',,si}=c6ab£'''s\ (7.21b) 
and all the other Poisson brackets vanish. Beware of the additional c factor appearing in the 



second bracket: our PN counting is different from that in [38|, and we have set S = cS'|js to 
keep close to the convention we use, see Eq. fll.ll) . One may as well employ a conserved-norm 
antisymmetric spin tensor defined (exactly) by: 



The dynamics is entirely contained in the Hamiltonian of the two-body system, expressed 
in terms of these canonical variables. Besides the recent result for the spin-orbit 3.5PN spin- 
orbit contributions jssf, the non-spin part of the Hamiltonian up to 2PN order can be found 



for instance in [63[, and the spin-orbit part up to 2.5PN order is given in [32[. As explained 
before, since we are working at 2PN relative order with respect to the leading order spin- 
orbit contribution, we will need to know the non-spin part of the dynamics only up to 2PN 
order. The time derivative of any quantity / is obtained as: 

f = {/.i/} + f. (7,23) 

where the second term accounts for a possible explicit time dependence of /, absent in our 
problem. Beware again that, because of the additional c factor in the Poisson brackets for 
the spin, the PN truncation of the Hamiltonian must be handled carefully. 

It is well known that, when comparing the results of the two formalisms, it is necessary 
to perform a contact transformation of the worldlines of the particles, i.e. a time-dependent 
shift of the worldlines that is not a global coordinate transformation. Thus, one must keep 
in mind that the positions and Xq differ at higher-order. We write this symbolically as 

y = Y(x,p,S). (7.24) 

The link between the harmonic velocity v and canonical momentum p is obtained by: 

v = V(x,p,S) = {Y(x,p,S),i/} . (7.25) 

As for the acceleration, we have: 

a = A(x,p,S) = {{Y(x,p,S),if},i/} . (7.26) 

On the other hand, we may use fl7.24p and (17.251) and the link between the spin variables 
(I7.30p to translate the harmonic-coordinates result (16. lip for the acceleration in terms of 
the ADM variables as: 

a = a(y,v,S) =a(Y(x,p,S),V(x,p,S),5(x,p,S)) , (7.27) 

and this expression must be identical to (17.260 . Thus, one is led to look for an extension of 
the contact transformation (I7.24p that realizes this identity. Focusing on the non-spin and 
spin-orbit parts, the structure of this contact transformation is the following: 

Yi = XI + i Y\-^™ + ^Y?™ + i Yf ™ + ^Y?™ + i Yf ™ + O ( \] . (7.28) 
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The lower-order already known expressions for this contact transformation can be found 



m 



63] and |32|. Notice that, as when using the Hamiltonian to compute time derivatives in 



the ADM setting, we need the non-spin part of this contact transformation up to 2PN only. 
To present these formulae, we use the canonical spin tensor variable f l7.22p to get the same 
index structure as in the other results of the present paper, and we adopt the convenient 
notation If a = Pa/'^^a shorten the expressions. We find: 
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77128 1 Tri + nil 
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1 ,^ 



--S2 7ri + {ni2n2)S2 + -(5'2ni27r2)n^2 



(7.29a) 
(7.29b) 

(7.29c) 



We need also the conversion rule between the harmonic spin tensor and the ADM canon- 
ical spin, say S = S(x.,p, S). The conversion rule from the conserved norm variable Spgg 
to S is given in Eqs. (6.4) and (6.5) of j32|, and the conversion rule between this conserved 
norm spin and Sfbb is given in Eq. (7.4) of 3l||. Finally, knowing the link between S 
and our spin tensor, Eq. (lAip in Appendix |Xl we get: 



FBB 



2PN ^ 
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^ij l_2_[j-r;j]^— fc 
^0,0 =^T^inSl 



(7.30) 

(7.31a) 
(7.31b) 
(7.31c) 



^0,1 = - {.ni2Ti2)n\2Si^A + 2(ni27r2)nf25'i'\2 + 47rf S'^^tt^ - 2(711211 2)1^ 2 S^i' 
-77r^5i 'k^+Ati'^S^ vr2 + (ni27r2) , 



71 



12 



"^0,2 ~'''12'-'l '''12 ^ '-''-'1 ' 

SA = — 8n^2'S'i ^^12 ~ Si 



(7.31d) 
(7.31e) 
(7.31f) 



Using the same structural hypothesis as in 32| (namely that in miYi terms such as 
rriiSi and 771282 are forbidden, where 71 is the total power of mass in the term), we build 
a putative spin contribution to the contact transformation with 112 unknown coefficients. 
Then, requiring the identity of the expressions (17.261) and fl7.27p above, we find a unique 
solution for these coefficients, among which 69 remain in the final result, which fixes uniquely 
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the contact transformation. Thus, we conclude that our result and the one obtained in the 
ADM formalism in Ref. H are equivalent^ 

We now give explicitly the 3.5PN spin-orbit extension of the contact transformation. 
For simplicity, we keep the same notation as we used for harmonic-coordinates results. In 
principle, we should write the result in terms of ADM variables, but since we are looking 
at the highest-order contribution here, only the leading order of variables makes sense and 
so all ADM and harmonic variables are equivalent (with v = p/m). We use the spin tensor 
variable, which at leading order is also equivalent to the ADM canonical spin tensor. The 
extension of the contact transformation reads: 



"^1 ^*3.5PN 
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We found a typographical error in the originally published version of the Hamiltonian (5) in Ref. [38|: the 
coefficient in front of the term _G_ ("i2Pi)P2 ICp, x P2)Si) should read —ts instead of — 
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VIII. CONCLUSION 

In this work, we computed the next-to-next-to-leading order spin-orbit contributions to 
the equations of motion of compact binaries. Those are of 3.5PN order for maximally 
spinning objects, thus improving our knowledge of the dynamics of such systems at this 
order. Our result was tested by checking the existence of a conserved energy, the manifest 
Lorentz invariance of the equations of motion, and the agreement of the test-mass limit 
with the motion of a spinless test particle around a Kerr black hole as well as that of a 
spinning test particle around a Schwarzschild black hole. We also recover and confirm the 
result obtained previously in Ref. [Isj using a reduced Hamiltonian method in ADM-type 
coordinates, extending the contact transformation that makes the link between the two 
formalisms. 

We leave for future work the study of the 3PN precession equation for the spins, the 
construction of the conserved quantities other than the energy (total angular momentum, 
linear momentum, center-of-mass integral), the center-of-mass reduction of the equations of 
motion as well as the further reduction of the dynamics to quasi-circular orbits. We shall 
also provide the spin-orbit contributions to the components of the near- zone metric itself and 
its value at the location of the particles. Most importantly, this work opens the way to the 
computation of the 3.5PN spin-orbit contributions to the energy flux emitted by the binary 
in the form of gravitational waves, and of the corresponding contribution in the phase of the 
GW signal, which should improve the templates used by current and future GW detectors. 



Appendix A: Link between different spin variables 

In this Appendix we provide the explicit conversion rule from the spin variable used in 
Paper I, wich we denote S'pgg, to our spin tensor variable S^^ . This rule allows one to readily 
translate the lower-order results of Paper I in terms of our variables for comparison. We 
adopt the notation (eaS) for e^^^a^ for any vector a, with the indices in this precise order: 



S { = \e^^^S{'' + \ S { + -, S i + 0(5) , (Al) 

FBB 2 IPN 2PN 

where 

S 1 = -\vle^^'S{' + \{ev,SM + ^e^^'^S{' , (A2a) 
IPN 4 2 2ri2 



2PN 16 4 ri2 
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12 
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(A2b) 



On the other hand the hnk between our spin tensor and the spin variable used in the ADM- 
Hamiltonian work 38| is provided in Eqs. f l7.30p - fl7.3ip . 
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